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Abstract 



O 

P i' 

g_) I We use t, 6, r Yukawa unification to constrain SUSY parameter space. We 

"-^ ■ find a narrow region survives for /i > (suggested by 6 ^ 57 and the 

anomalous magnetic moment of the muon) with Aq ~ —1.9 m-ie, "^lo ~ 

|X| ■ 1-4 mi6, mi6 ~ 1200-3000 GeV and/x,Mi/2 ~ 100-500 GeV . Demanding 

H ! Yukawa unification thus makes definite predictions for Higgs and sparticle 

masses. 



Minimal supersymmetric [SUSY] SO (10) grand unified theories [GUTs] have many 
profound features ^: all fermions in one family sit in one 16 dimensional spinor rep- 
resentation; the two Higgs doublets of the minimal supersymmetric standard model sit 
in one 10 dimensional fundamental representation, and gauge coupling unification at a 
GUT scale Mg ~ 3 x 10^^ GeV fits well with the low energy data g |. In addition 
in the simplest version of SO(IO) the third generation Yukawa couplings are given by 
a single term in the superpotential PF = A 16 10 16 resulting in Yukawa unification 
A* = Aft = At- = Aj,^ = A and a prediction for Mt with large tan/? ~ 50 [Q. [] This beau- 
tiful result is however marred by potentially large weak scale threshold corrections [0, H 



"^ -" ^ 9 _L A^X+ , A^x" , A^logN 



TTibiMz) = XbiMz) -^ cosjS (1 + Ami + ^^b + ^^b + ^^, 

For /i > the gluino term is positive and in most regions of SUSY parameter space it 
is the dominant contribution to Arrib . Reasonable fits prefer Arub < 0; hence Yukawa 
unification is easy to satisfy with /x < 0. 

The decay b —>■ s'j and the muon anomalous magnetic moment also get significant 
corrections proportional to tan/3[^. The SUSY contribution to b ^ s'j comes from one 
loop diagrams similar to those contributing to the bottom mass. The chargino term 
typically dominates and has opposite sign to the SM and charged Higgs contributions, 
thus reducing the branching ratio for /i > 0. This is necessary to fit the data since the 
SM contribution is somewhat too big. fi < would on the other hand constructively 
add to the branching ratio and is problematic. In addition, the recent measurement of 
the anomalous magnetic moment of the muon a^^^ = {q ~ 2)/2 = 43 (16) x 10~^° also 
favors /i > [0|. Thus it is important to confirm that Yukawa unification can work 
consistently with /i > . 

In this paper we assume exact Yukawa unification and search, using a x^ analysis, 
for regions of SUSY parameter space with fi > providing good fits to the low energy 
data. We show that Yukawa unification dramatically constrains the Higgs and SUSY 
spectra. These results are sensitive to the SUSY breaking mechanism. 

It is much easier to obtain EWSB with large tan/3 when the Higgs up/down masses 
are split {'m'jj^ < mfj ) [^ . In our analysis we consider two particular schemes we refer to 
as universal and D term splitting. In the first case the third generation squark and slepton 
soft masses are given by the universal mass parameter rriiQ , and only Higgs masses are 
split: m^H^ HA = rnlQ (1 =F Arriff). In the second case we assume D term splitting, i.e. 
that the D term for U{l)x is non-zero, where U{l)x is obtained in the decomposition 
of 5*0(10) — > SU{5) X U{l)x- In this second case, we have rn^Hu,Ha) ~ "^lo "F 2-Dx, 
^(Q, u, e) = ^w + Dx, mfj^ L) = ^l& 



mfn y^ e) ~ ""^16 + ^x, ^j T) = "^16 ~ 3L'x- The universal case does not at first sight 



""^Note, GUT scale threshold corrections to this Yukawa unification boundary condition are naturally 
small ( < 1% ), since they only come at one loop from the SO(IO) gauge sector and the third generation 
- Higgs Yukawa coupling B. This is in contrast to GUT scale threshold corrections to gauge coupling 
unification which may be significant, coming from doublet/triplet splitting in the Higgs sector and, even 
more importantly, the SO(IO) breaking sector which typically has many degrees of freedom. The data 
requires eg - aoiMc) 4/o g. 
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Figure 1: x^ contours for mie = 1500 GeV (Left) and mig = 2000 GeV (Right). The 
shaded region is excluded by the chargino mass hmit m^+ > 103 GeV. 



appear to be similarly well motivated. It is quite clear however that in any SUSY model 
the Higgs bosons are very special. R parity is used to distinguish Higgs superfields 
from quarks and leptons. In addition, a supersymmetric mass term /i with value of 
order the weak scale is needed for an acceptable low energy phenomenology. Since /i is 
naturally of order Mc, one needs some symmetry argument why it is suppressed. Of 
course, if the Higgs are special, then it is reasonable to assume splitting of Higgs, while 
maintaining universal squark and slepton masses. This may be achieved by GUT scale 
threshold corrections to the soft SUSY breaking scalar masses |jTO[ . Here we present the 
most compelling mechanism |]ll|. In 5*0(10), neutrinos necessarily have a Yukawa term 
coupling active neutrinos to the "sterile" neutrinos present in the 16. In fact for v^- we 
have Xy^ Vr L Hu with Xy^ = Xt = Xb = Xr = X. In order to obtain a tau neutrino with 
mass m^^ ~ 0.05 eV (consistent with atmospheric neutrino oscillations), the "sterile" 
Ut- must obtain a Majorana mass Mp^ > 10^^ GeV. Moreover, since neutrinos couple to 
Hu (and not to Hd) with a fairly large Yukawa coupling (of order 0.7), they naturally 
distinguish the two Higgs multiplets. With A = 0.7 and Mp^ = 10^^ GeV, we obtain a 
significant GUT scale threshold correction with Am^ ^ 10%, remarkably close to the 
value needed to fit the data. At the same time, we obtain a small threshold correction 
to Yukawa unification ^ 2.5% (for more details see 0]). 

Our analysis is a top-down approach with 11 input parameters, defined at Mg, 
varied to minimize a x^ function composed of 9 low energy observables. The 11 in- 
put parameters are: Mq, aciMc), €3; the Yukawa coupling A, and the 7 soft SUSY 



breaking parameters fi, M1/2, Aq, tan/3, m 



.16, mJQ, Amjj {Dx] 



for universal (D term) 
case. We use two (one) loop renormalization group [RG] running for dimensionless (di- 
mensionful) parameters from Mq to Mz and complete one loop threshold corrections 
at Mz |§]. We require electroweak symmetry breaking using an improved Higgs po- 
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Figure 2: Contours of constant mh{'mh)[GeV] (Left) and Arrify in 
2000 GeV . 
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tential, including mf and m^ corrections in an effective 2 Higgs doublet model below 
Mstop [|12]. The x^ function includes the 9 observables; 6 precision electroweak data 
Mz, Mw, pNEW and the 3 fermion masses Mtop, mf,{mb), Mr. The 



a^M, G^, OisiMz 



experimental values used for the low energy observables are given in the table. 

Fig. 1 shows the constant x^ contours for rriiQ = 1500 and 2000 GeV in the case of 
universal squark and slepton masses. We find acceptable fits (x^ < 3) for Aq ~ — 1.9 tjiiq, 
rriio ~ 1-4 rriiQ and itliq > 1.2 TeV. The best fit is for itliq > 2000 GeV with x^ < 1- 
Note, electroweak symmetry breaking in this region of parameter space requires split- 
ting Higgs up/down masses, Amj^ ~ 0(13%). This range of soft SUSY parameters is 
consistent with solution (B) of Olechowski and Pokorski ||T0|. In the table we present 
the input parameters, the fits and the predicted Higgs and SUSY spectra for two repre- 
sentative points with universal squark and slepton masses and the best fit value for D 
term splitting. We have not presented the contour plots for D term splitting since as 
can be seen from the best fit point in the table, the bottom quark mass is poorly fit in 
this case and x^ > 5. Recall, since we have 11 input parameters and only 9 observables, 
we consider such poor fits unacceptable. 

Fig. 2 gives the constant mbirrib) and Anib contours for rriiQ = 2000 GeV . We 
see that the best fits, near its central value, are found with Arub < —2%. Why does 
Yukawa unification only work in this narrow region of SUSY parameter space? The log 
corrections Am^^ ~ 4 — 6% (total contribution from gluino, neutralino, chargino and 
electroweak loops) are positive and they must be cancelled in order to obtain Anib < —2 
%. The leading mass insertion corrections proportional to tan/3 are approximately given 
by0 

^ tan/3 and Am^ ^ * — 2" tan/S. 
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Figure 3: Contours of constant vtih [GeV] (solid lines) with x^ contours from Fig. 
(dotted lines) for mig = 1500 GeV (Left) and mig = 2000 GeV (Right). 



They can naturally be as large as 40%. The chargino contribution is typically opposite 
in sign to the gluino, since At runs to an infrared fixed point oc — Mi/2(see for example, 
Carena et al. [^). Hence in order to cancel the positive contribution of both the log 
and gluino contributions, a large negative chargino contribution is needed. This can 



be accomplished for —At > nig and m^- 



« ^b,- 



The first condition can be satisfied 



for Aq large and negative, which helps pull At away from its infrared fixed point. The 
second condition is also aided by large At. However in order to obtain a large enough 
splitting between m^^ and m^ , large values of mig are needed. Note, that for universal 
scalar masses, the lightest stop is typically lighter than the sbottom. We typically find 
mi_^ ~ 3 rrit^. On the other hand, D term splitting with Dx > gives m^^ < mj^. As 
a result in the case of universal boundary conditions excellent fits are obtained for top, 
bottom and tau masses; while for D term splitting the best fits give mb{rnh) > 4.59 GeV . 
Finally in Fig. 3 we show the constant light Higgs mass contours for rriiQ = 1500 and 
2000 GeV (solid lines) with the constant x^ contours overlayed (dotted lines). Yukawa 
unification for x^ ^ 1 clearly prefers light Higgs mass in a narrow range, 112 - 118 GeV . 
In this region the CP odd, the heavy CP even Higgs and the charged Higgs bosons are also 
quite light (see fit 2 in the table). [| In addition we find the mass of ti ~ (150 — 250) GeV, 
b^ ^ (450 - 650) GeV, n ~ (200 - 500) GeV, g ~ (600 - 1200) GeV, x+ ~ (100 - 250) 
GeV, and x^ ~ (80 — 170) GeV. All first and second generation squarks and sleptons 

^ It would be interesting to see how sensitive our results, for Higgs masses, may be to alternative 
electroweak symmetry breaking approximations. In this paper we have used the effective 2 Higgs doublet 
analysis of |I3| with an estimated 3 GeV uncertainty in Higgs masses. This approximation may be 
particularly well suited to the light Higgs spectrum we obtain in our analysis. The alternative scheme, 
in which the Higgs tadpoles are evaluated at a scale of order M^top Iql is however more frequently used 
in the literature. 
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Table 1: Three representative points of the fits. 

We fit the central values: Mz = 91.188, Mw = 80.419, G^ x 10^ = 1.1664, a 
137.04, Mj. = 1.7770 with 0.1% numerical uncertainties; and the following with the 
experimental uncertainty in parentheses: as{Mz) = 0.1180 (0.0020), pnew x 10^ = 
-0.200 (1.1), Mt = 174.3 (5.1), mb{mb) = 4.20 (0.20). The neutral Higgs masses 
h, H, Aq are pole masses; while all other sparticle masses are running masses. 
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Input parameters 








aa' 


24.46 


24.66 


24.73 


Mg X 10-16 


3.36 


3.07 


3.13 


es 


-0.042 


-.0397 


-0.046 


A 


0.70 


0.67 


0.80 


miQ 


1500 


2000 


2000 


mio 


2027 


2706 


2400 


Amjj 


0.13 


0.13 


0.07 


Ml/2 


250 


350 


350 


P 


150 


200 


115 


tan/3 


51.2 


50.5 


54.3 


^0 


-2748 


-3748 


-731 


X^ observables 








Mz 


91.13 


91.14 


91.15 


Mw 


80.45 


80.45 


80.44 


G^ X 10^ 


1.166 


1.166 


1.166 


O^EM 


137.0 


137.0 


137.0 


as{Mz) 


0.1175 


0.1176 


0.1161 


Pnew X 10^ 


0.696 


0.460 


0.035 


Mt 


175.5 


174.6 


177.9 


mb{mb) 


4.28 


4.27 


4.59 


Mr 


1.777 


1.777 


1.777 


TOTAL x^ 


1.50 


0.87 


5.42 


ho 


116 


116 


115 


Ho 


120 


121 


117 


Ao 


110 


110 


110 


H+ 


148 


148 


146 


X? 


86 


130 


86 


xl 


135 


190 


126 


xt 


123 


178 


105 


9 


661 


913 


902 


h 


135 


222 


1020 


h 


433 


588 


879 


n 


288 


420 


1173 


al^^' X W 


9.7 


5.5 


6.1 



have mass of order rriie- The hght stop and chargino may be visible at the Tevatron. 
With this spectrum we expect ti ^ x^ b with x^ ^ x^ I u to be dominant. Lastly Xi 
is the LSP and possibly a good dark matter candidate [0. 

The region of SUSY parameter space preferred by Yukawa unification may be consis- 
tent with a supergravity mechanism for SUSY breaking at Mpi with RG running from 



Mpi to Mq (see for example Murayama et al. |T0|). It however cannot be obtained with 
gauge mediated or gaugino mediated SUSY breaking mechanisms where Aq = at zeroth 
order. It may also be obtained in anomally mediated schemes but in this case one still 
has to worry about slepton masses squared and also the fact that in this case, since the 
gluino and chargino masses have opposite sign, it is difficult to fit both 6 — i> 57 and a^. 



In a future paper [|TT| we present the sparticle spectrum in more detail and conse- 
quences for Tevatron searches. We discuss the sensitivity of our results to small GUT 
scale threshold corrections to Yukawa unification with both universal and D term Higgs 
up/down splitting. We also check the robustness of the Higgs spectrum by artificially 
adjusting the CP odd Higgs mass using a penalty in x^. We find that x^ increases by 
at most 40% for any m^o less than ^ 350 GeV. The light Higgs mass rrih is rather 
insensitive to the value of rriAO', whereas ■mH,TnH+ are linearly dependent on m^o. We 
also consider constraints resulting from the processes b — »• sj, Bg -^ fi~^ fi~ , a^^'^ and 
the proton lifetime in a semi model independent way. We shall only make a few short 
comments here. In order to fit 6 — »• 57 we find that the coefficient C^^^^^'^ ~ —C^^ (see 
for example, Eqn. 9 in Ref. |T^) with the chargino term dominating by a factor of order 



5 over all other contributions. This is due to the light stop ti. In fact, 6 — i> 57 is more 
sensitive to mi^ than mb{mb). Fitting the central value B{b -^ 57) = 2.96 x 10~^ [IT3 



requires a heavier ti with {mi^)MiN ~ 500 GeV; significantly larger than the range which 
provides the best fits to m;,. We now find mb{mb)MiN ~ 4.3. Moreover no other sparticle 
masses are affected. The process Bg -^ fi^ fi^ provides a lower bound on rrij^o > 200 
GeV (see recent work of [|16|). However this has only a minor impact on x^ as discussed 



above. We recall that proton decay experiments prefer values of miQ > 2000 GeV and 
rriiQ » Ml/2 (see ref. fl^). This is in accord with the range of SUSY parameters found 
consistent with third generation Yukawa unification. There is however one experimental 
result which is not consistent with either Yukawa unification or proton decay and that 
is the anomalous magnetic moment of the muon. Large values of rriiQ > 1200 GeV lead 
to very small values for a^^^ < 16 x 10~^°. Hence a necessary condition for Yukawa 
unification is that forthcoming BNL data P| and/or a reanalysis of the strong interaction 
contributions to a^^ will significantly decrease the discrepancy between the data and 
the standard model value of a^. 

In summary, most of the results of our analysis including only third generation 
fermions remain intact when incorporating fiavor mixing. The light Higgs mass and most 
sparticle masses receive only small corrections. The lightest stop mass increases, due to 
6 — > 57. Nevertheless there is still a significant ti — ^2 splitting and m^-^ << m^ . The 
A^, H, H^ masses are necessarily larger in order to be consistent with B^ —>■ /x"*" /x~ [ITT 



■^We thank K.S. Babu and C. Kolda for discussfons. 



which suggests that this process should be observed soon; possibly at Run II of the 
Tevatron. Finally, the central value for a^^'^ must significantly decrease. The "smok- 
ing guns" of SO(IO) Yukawa unification, presented in this letter, should be observable 
at Run III of the Tevatron or at LHC Also, in less than a year we should have more 
information on a^,^'^ . 

fj. 

In previous works Yukawa unification with /i > was not possible.^ Pierce et 
al. [§] assume Am^ = and, as a result, they are not able to enter the region of SUSY 



parameter space consistent with both EWSB and Yukawa unification. Baer et al. [jl8 
also cannot obtain Yukawa unification with /i > . This is because they use D term 
splitting for Higgs up/down which as discussed typically leads to sbottom lighter than 
stop. 

While completing this article, the paper by Baer and Ferrandis [|19| appeared which 



confirmed our results [^ on the existence of a preferred region of SUSY parameter space 
consistent with Yukawa unification and /i > . Their results however require significant 
GUT threshold corrections to Aj = A^ = A,- of order 8 - 15% which helps them obtain 
mi^ < mi . They also claim that better fits are obtained with D term splitting than 
with the universal splitting case. We believe the latter is only true because the authors 
do not allow their SUSY parameters, in particular rriiQ and Aq, to explore the region of 
parameter space discussed in |^ and this paper. Q 
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